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^ : Abstract 
cn 

We estimate terrestrial antineutrino and neutrino fluxes according to dif- 

CN . ferent models of Earth composition. We find large variations, corresponding 

I to uncertainties on the estimated U, Th and K abundances in the mantle. 

' Information on the mantle composition can be derived from antineutrino flux 



o 



measurements after subtracting the crust contribution. This requires a good 
description of the crust composition in the region of the detector site. Mea- 



. surements of terrestrial antineutrinos will provide a direct insight on the main 

'. sources of Earth's heat flow. 

V 

■ I. INTRODUCTION 

Earth emits a tiny heat flux with an average value $h = 80mW/m'^ , deflnitely smaller 
than the radiation coming from the Sun, Kq = lAkW/m?, larger however than the energy 
^ \ deposited by cosmic rays, $c — l^~^W/m?. When integrated over the Earth surface, the 
c3 ■ tiny flux translates into a huge heat flow, H^^ ~ AQTW , the equivalent of ten thousand 
nuclear power plants [1]. 

We would like to recall to the particle physics community that the sources of Earth energy 
flow are not understood quantitatively and that measurements of (anti)neutrinos from the 
Earth in the next few years should be capable of determining the radiogenic contribution. 

A comparison between the Sun and Earth energy inventories may be useful for illustrating 
the differences in the two cases. Clearly, a heat flow H can be sustained for a time t provided 
that an energy source of at least U = Ht is available. 

For the Sun, U = Hq ~ 5 ■ 10^^ J and clearly neither gravitation {Ug — GMq/Rq = 
4 ■ 10^^ J ) nor chemical reactions {Uch — 0.1 eV ■ Nq = 2 ■ 10^'' J, where is the number of 
nucleons) are enough, and only nuclear energy {Unuc — 1 MeV Nq = 2 ■ 10^^ J) can sustain 
the solar luminosity over the solar age, as beautifully demonstrated Gallium experiments in 
the last decade [2]. On the other hand for the Earth one has Ug — 4:- 10^^ J, Uch — 6 ■ 10^^ J 
and Unuc — 6 ■ 10^° J (assuming some that some 10~^ of Earth mass consists of radioactive 
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nuclei), so that each of the previous mechanisms in principle can account for = 5 ■ 10^'' J. 
In order to understand the energetics of the Earth one has to clarify the roles of the different 
energy sources, their locations and when they have been at work. At the end of a review on 
the Earth energy sources J. Verhoogen [3] summarized the situation with the following words: 
"What emerges from this morass of fragmentary and uncertain data is that radioactivity by 
itself could plausibly account for at least 60 percent, if not 100 percent, of the earth's heat 
output. If one adds the greater rate of radiogenic heat production in the past, . . . possible 
release of gravitational energy (original heat, separation of core, separation of inner core, 
tidal friction,... meteoritic impact ...), the total supply of energy may seem embarrassingly 
large. . . . Most, if not all of the figures mentioned above are uncertain by a factor of at least 
2, so that disentangling contributions from the several sources is not an easy problem. ". 

In this respect a determination of the radiogenic contribution is most important. Ra- 
diogenic heat arises mainly^ from the decay (chains) of ^^^t/, ^^^T/i and '^^K. All these 
elements produce heat together with antineutrinos, with well fixed ratios heat /neutrinos. A 
measurement of the antineutrino flux, and possibly of the spectrum, would provide a direct 
information on the amount and composition of radioactive material inside Earth and thus 
would determine the radiogenic contribution to the heat flow. 

On the other hand, until recently the neutrino fate could not be predicted reliably, as 
testified by the thirty years old solar neutrino puzzle [4] . The disagreement between theory 
and observation by factors of order two suggested that (anti) neutrino survival probabilities 
were essentially known within factors of two. Thus observation of terrestrial (anti) neutrinos 
could not be useful for improving our knowledge of Earth radioactivity. The situation 
has dramatically changed since the SNO results [6], which clearly prove that a fraction of 
electron neutrinos change their flavor during the trip form Sun to Earth. When combined 
with the results of other solar and terrestrial neutrino experiments, the picture is converging 
towards the so called large mixing angle (LMA) oscillation solution. In other words, now 
we can predict reliably the fate of terrestrial neutrinos and antineutrinos, in their trip from 
production site to detectors. 

Last but not least, the experimental techniques for detection of MeV antineutrinos have 
enormously improved in the last few years. As testified by the development of Kamland 
[7] and Borexino [8], it is now possible to build kiloton size detectors, with extremely low 
background. 

The argument of geo-neutrinos was introduced by Eder [9] in the sixties, it was reviewed 
extensively by Krauss, Glashow and Schramm [10] in the eighties and it has been considered 
more recently in [11,12]. Now it is the right time for neutrino physics to contribute in 
reconstructing the thermal history of the Earth. 



For simplicity we neglect ^^°U and Rb which provide smaller contributions. 
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II. ENERGY SOURCES AND NEUTRINO LUMINOSITIES 



The heat production rates per unit mass of natural U, Th and K are given by^: 

e(C/) = 0.95 • 10"^ W/kg (1) 
e(T/i) = 0.27 • 10-^ W/kg 
e(ir) = 0.36 • 10-^ W/kg 

This is sufficient to determine the Earth radiogenic heat production rate H in terms of 
the mass of each element. When heat production is expressed in TW and masses in units 
of 10^^ kg one has: 

H = 9.5 M{U) + 2.7 M{Th) + 3.6 • 10"^ M{K) (2) 

It is convenient to write this equation in terms of the uranium mass M{U) and of the mass 
ratios of the other elements to [/, as these latter quantities are more regularly distributed 
in terrestrial and meteoritic samples: 

H = 9.5 M{U) [1 + 0.28 Th/U + 3.8 • 10"^ K/U] (3) 

The specific neutrino production rate (neutrinos per unit mass and time) of each element 
e^,, is immediately derived from the isotopic abundance, decay time and the number of 
neutrinos emitted in each decay, see Table I. (Anti) Neutrinos luminosities are immediately 
derived in terms of the mass of each element and the appropriate Cp/j,. Measuring Lpj^, in 
units of 10^^ particles per second and masses in units of 10^'' kg one has: 

Lp = 7.4 M{U) + 1.6 M{Th) + 2.7 • 10"^ M{K) (4) 
= 3.3- lO-^M(X) 

(5) 

We have thus the basic equations for determining radiogenic heat production and neu- 
trino flows from models of the Earth composition. 

A. A naive chondritic earth 

The simplest model assumes that the global composition of the Earth is similar to that 
of the oldest meteorites, the carbonaceous chondrites (CI). 

The typical values of CI [5] are Th/U = 3.8 , K/U = 7-10^ and U / Si = 7.3 ■ 10"^ 
[1]. Silicon represents about 15% of the Earth mass, Af© = 5.97 ■ 10^^ kg. If these elements 



^The marked difference of e{K) corresponds to the fact that the natural abundance of '^^K is 
1.2 • 10-'', i.e. ei'^^K) = 0.3 • 10^ W/kg. 
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have not been lost in the Earth formation process, one obtains M(U) — 0.653 • lO-*^^ kg, 
M{Th) = 2.48 ■ 10^^ kg and M{K) = 4.57 ■ lO^^ kg. 

The contribution to heat flows and neutrino luminosities are reported in the flrst column 
of Table II. Radiogenic production in the chondritic model easily accounts for 75% of the 
observed heat flow, and it could easily saturate it when uncertainties are included. Uranium 
and thorium provide comparable contributions, each a factor of two below that of potassium. 
Concerning antineutrinos, potassium dominates by an order of magnitude at least, as a 
consequence of the more favourable neutrino/energy ratio. 

B. The Bulk Silicate Earth model 

Uranium, thorium and potassium are lithophile elements, so they should accumulate in 
the Earth crust. Actually, upon averaging data over the huge differences between continental 
and oceanic components, it is found that Earth crust contains some 3/4 of the uranium 
predicted for the whole Earth by the chondritic model [1]. Within large variations, Th/U 
is consistent with the chondritic prediction. On the other hand, the crust looks depleted in 
potassium, the typical ratio being K/U — 10, 000, a factor 7 below that of CI. 

Observational data on the mantle, which are anyhow limited to the upper part, suggest 
that uranium and potassium are globally more abundant than the CI prediction, Th/U is 
consistent with the chondritic value and the potassium depletion is conflrmed. No observa- 
tional data are available on the core, which should consist of siderophile elements without 
signiflcant amount of U, Th or K. 

Actually, when deriving Earth composition from mctcoritic data, one has to take into 
account the volatilization of a signiflcant fraction (some 17%) [15] of the total Si02, so that 
a larger amount of metoritic material is needed for Earth formation. 

The origin of potassium depletion, also observed in the Moon, Venus and Martian me- 
teorites, is somehow uncertain. Elements of the atomic weight of potassium cannot be lost 
from the terrestrial planets, even at elevated temperatures, once these bodies have reached 
their present size [16]. The most reasonable explanations seems that this clement was de- 
pleted in the precursor planctcsimals from which the inner planets accumulated^. 

All this brings us to the Bulk Silicate Earth (BSE) model, which provides a description 
of geological evidence coherent with geochemical information. It describes the primordial 
mantle, prior to crust separation. The estimated uranium mass is: 

M(U) =0.8-10^^ kg , (6) 

within some 20% [17], the ratio Th/U is close to the chondritic value and K/U = 10,000. 
The present crust and mantle should contain respectively about one half of each element. 



^It has been suggested that potassium behaves as a metal at high pressure, and thus it can be 

buried in the planetary cores. This hypothesis could work for Earth, and it provides a suitably 
placed energy source for sustaining the terrestrial magnetic field, see [1]. However it does not 
explain potassium depletion in Mars, where the central pressure, only 400 kbars, is insufficient for 
potassium to enter a Martian core. 
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In this BSE model the (present) radiogenic production, mainly from uranium and tho- 
rium, accounts for about one half of the total heat flow. The antineutrino luminosities from 
uranium and thorium are rescaled by a factor 1.3 whereas potassium, although reduced by 
a factor of 5, is still the principal antineutrino source. 

C. A fully radiogenic model 

At the other extreme, one can conceive a model where heat production is fully radiogenic, 
with K/U fixed at the terrestrial value and Th/U at the chondritic value, which seems con- 
sistent with terrestrial observations. All the abundances are rescaled so as to provide the full 
40 TW heat flow (last column of table II). All particle production rates are correspondingly 
re-scaled by a factor of two with respect to the predictions of the BSE model. 



In summary, the discussion of these somehow extreme models shows that particle lumi- 
nosities are uncertain by a factor of order two, the relative contributions to heat production 
are strongly model dependent whereas potassium is anyhow the principal neutrino and anti- 
neutrino source. 



III. FROM LUMINOSITY TO FLUX AND SIGNAL 

An (anti) neutrino detector near the Earth surface {R = i?®) is sensitive to the flux 
impinging onto it from any direction: 

where the integral is taken over the Earth volume and A is the number of particles produced 
per unit volume and time^. The flux depends on the geometrical distribution of the sources, 
and we we can write: 

where G is a geometrical factor of order unity. One has: $/(10^ cm~^ s"^) = 
0.2GL/{10^^s-^). 

In order to estimate fluxes, one needs to know the distribution of radioactive elements 
in the Earth interior, and not only their total abundances. Concerning the crust, these 
elements are mainly concentrated in the continental part. Recent estimates of the uranium 
average mass abundance in the continental crust (CC) are near acc{U) = 1.7 ppm [13]. The 



^This is different from the flux normal to the Earth surface, which in the case of spherical symmetry 
is given by $_l = J (frA{r). 
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abundance in the oceanic crust (OC) is an order of magnitude smaller, aoc — 0-1 ppm. If one 
also considers that OC is much thinner than CC [Mcc = 2.3-10^^ kg and Mqc = 0.6-10^^ kg) 
one concludes that the contribution of the oceanic crust to neutrino production is definitely 
smaller. In this way we estimate the total uranium mass in the crust: 

Mc{U) ^ OA ■ 10^^ kg . (9) 

For the sake of a first estimate, we shall assume that Mc{U) is uniformly distributed 
over the Earth surface within a layer of thickness h — 30 km. We associate to the mantle 
an uranium mass Mjn{U) = M{U) — Mc{U), where M{U) was estimated in the previous 
section for each model and Mc{U) is given from eq.(9). Furthermore, we shall assume 
uniform distribution within the mantle. We also assume a uniform distribution of Th/U 
and K/U. Within these approximations the geometrical factors G are easily calculated. 

For a spherical shell with radii ri = XiR^ and r2 = 2^2-^0 and uniform distribution one 
has: 

«-2M^{i(^-^')'"^-^'-i(^-^Slni±f + .4 (10) 

For the crust (r2 = -R® and r2 — ri = /i ~ 30 km) and for the mantle (r2 — -R® and 
ri ~ -Re/2) one has: 

G'e~^[l + M2i?eA)] = 3.5 ; (7^ = 1.6 (11) 

This allows the calculation of the fluxes for each Earth model reported in Table III. 

One remarks that fluxes arc of the order of magnitude of the solar Boron neutrino 
flux, potassium (anti)neutrinos are the dominant component in any model. The various 
models yield signiflcantly different predictions. Contributions from crust and mantle look 
comparable. 

Indeed in order to reach these results we used several approximations, which is worth 
discussing, also in view of obtaining more precise estimates. 

i) Global effects. 

There are signiflcant uncertainties on the value of acc{U). For the lowest value found in the 
literature, acc{U) = 0.91 ppm [14], the uranium mass in the crust is halved with respect 
to the reference value in eq.(9). In addition, there are indications that the upper part 
of the mantle is impoverished in the content of radioactive elements. The effect of these 
uncertainties are shown in Table IV, where ^p{U) is calculated for a flxed (BSE) value of 
the crust +mantle uranium mass. A change of the [/-abundance in the crust by a factor 
of two corresponds to a 15% change in the flux. If uranium is completely removed from 
the upper mantle the flux is reduced by 3%. All these effects are thus smaller with respect 
to the uncertainty on the total uranium mass, which is reflected in changes of fluxes by 
factor two when going through the extreme chondritic and fully radiogenic models. Similar 
considerations hold for thorium and potassium 

ii) Regional effects. 

The fluxes reported in Table III are averaged values and one has to remind that Earth crust 
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is significantly variable in thickness and composition. Thus one has to expect significant 
variations of the actual fluxes, depending on the detector location. 

As a few significant examples, the Gran Sasso laboratory sits on a thick continental 
crust, whereas Kamioka is on an island arc in between the Eurasian plate and the Pacific 
plate. Following the approach of ref. [12] we have estimated $p/,^ at the two sites, in a model 
which distinguishes between CC and OC and includes a variable crust thickness. The model 
is based over a global crustal map at 5° x 5", and assumes imiform distribution within the 
mantle. Uranium mass in the crust and in total Earth are kept fixed at the values of eqs. (6) 
and (9). Table V shows that fiuxes at specific sites can differ significantly from the average 
value and thus a crust map is really needed for a precise fiux estimate at the detector size. 
On the other hand, the difference between Gran Sasso and Kamioka is at the level of 15%. 
Local effects. 

The flux from the crust within a distance d from the detector is easily estimated using planar 
geometry [d « i?©): 

$(< d) = ^[d arctan(/i/d) + {h/2)ln{l + (^/h^)\ (12) 

where h is the local crust thickness and A is the local activity. This has to be compared 
with the total fiux from the crust: 

= GM- (13) 

where A and h are the mean crustal activity and thickness. 

A significant quantity is the relative contribution R = $(< d)/^^ ■ By suitable expan- 
sions of eq. (12) one immediately derives the contribution of the nearby rocks (say d = h/3) 
and of the local area, i.e. up to a distance d « i?©: 

Rnea = 0.07Ah/{Ah) (14) 

Rioc = (1/7)[1 + ln{d/h)]Ah/{Ah) . (15) 

(16) 

H A — A and h — h — 30 km one finds that the rocks within 10 km contribute 7% 
with respect to the total from the crust. Within 100 km one has about 30% of the crust 
contribution and about 15% of the total flow, from crust and below. A more detailed 
geochemical information of this area is thus needed if one aims at a few percent accuracy 

In summary, detection of terrestrial (anti)neutrinos is particularly important for deter- 
mining the amount of radioactive material inside the mantle, where information on the 
chemical composition is most uncertain, the lower part being completely unaccessible to 
observation. A suitable approach would thus consist of i) measuring the (anti) neutrino 
flux; ii) subtracting the component originated in the crust, which has been mapped with 
geological methods, so as to determine the corresponding abundances in the mantle. The 
previous calculations show that crust and mantle provide comparable fiuxes, so that the 
subtraction procedure is possible. The crust description however has to be more detailed in 
the proximity of the detector. 

A detailed discussion of the (anti)neutrino signal is beyond the aim of this letter and we 
would like to remark just a few relevant points. 
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i) Due to the different antineutrino energy end-points [9] {Ej^ax— 3.26, 2.25 and 1.31 
MeV for U, Th and '^^K respectively) it is possible at least in principle to separate the 
contributions to Earth radioactivity. 

ii) Antineutrinos from U and Th can be detected and separated by means of 

P + p ^e+ + n- 1.804 MeT/ , (17) 

whereas different detection schemes are necessary for K antineutrinos, which are below the 
energy threshold for (17). The monochromatic {E = 1.513 Mel^ ) neutrinos from ^^K are 
essentially obscured by the dominant solar flux {^(pep) — 1.4-10^ cm~^ at E — 1.44 MeV 
and ^(CNO) ~ 10^ cm-^ s'^ MeV . 

iii) So far we did not consider the effect of neutrino oscillations. If LMA is the correct 
solution, with Am^ ~ 5.5 • 10'^ eV^ and sm^2e ~ 0.83 [18], the oscillation length L = 
4:71 E/Am'^ of 1 MeV antineutrino is around 45 km and the (distance averaged) survival 
probability of electron antineutrinos Pge = 1 — l/2sin^2^^ ~ 0.58 The present uncertainty 
on sin^ 29 (about 20%) translates into a 15% uncertainty on the fluxes. 

iv) Events from U and Th antineutrinos in a organic scintillator detector have been 
estimated in the range (20-100)/kton-year [11,12], so that a flux measurement with a 10% 
accuracy should be feasible in a few years. The main background source [19] is antineutrinos 
from nuclear power plants (see last row of Table III), which depends on the detector location. 

IV. CONCLUDING REMARKS 

We have estimated terrestrial antineutrino and neutrino fluxes according to different 
models of Earth composition. We flnd large variations, corresponding to uncertainties on 
the estimated U, Th and K abundances in the mantle. Information on the mantle composi- 
tion can thus be derived from (anti) neutrino flux measurements after subtracting the crust 
contribution. This requires a good description of the crust composition in the region of the 
detector site and in return it will provide direct insight on the main sources of Earth's heat 
flow. 

Just a few years after the celebrated slow neutron studies of the Rome group, Bruno 
Pontecorvo developed the neutron well log [20] , an instrument which is still used in geology 
for the search and analysis of hydrogen containing substances (water and hydrocarbons). 
Possibly it is now the time for applying to different disciplines what we have learnt so 
far on neutrinos. In fact, there are several attempts in this direction, see e.g. [21] and 
references therein. The determination of the radiogenic component of the terrestrial heat is 
an important and so far unanswered question. It looks to us as the flrst fruit which we can 
get from neutrinos, and Kamland will catch the flrstlings very soon. 
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Addendum 

After this paper was submitted, the first results of KamLAND have become available [22]. 
From an exposure of 1.39 ■ 10^^ protons-year, 9 geo-neutrino events are reported. For the 
best fit survival probability ^(t'e ^e) = 0.55 [22], from the average fiuxes calculated in 
Table III one predicts 2.6, 3.1 and 5.1 events respectively for the chondritic, BSE and fully- 
radiogenic models. Predictions from fiuxes calculated specifically for the Kamioka site are 
3.5, 4.0 and 6.0 events respectively. They are all consistent with the experimental value, 
within its statistical fiuctuation of ±5.7 events [23]. 
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TABLES 



TABLE I. Main radiogenic sources. We report the Q-values, the half hves (ri/2), the 
maximal energies {Emax) and (anti)neutrino production rates (ep/,^) per unit mass for natural 
isotopic abundances. Neutrinos from electron capture are monochromatic. 



decay 


Q 


n/2 


Emax 








[MeV] 


[109 y] 


[MeV] 


[kg-' 




238 ^206 p5 ^ 84^g + 6e + 6i> 


51.7 


4.47 


3.26 


7.41 • 


10^ 


232 j.^ ^208 + g4^g + 4e + 4Z^ 


42.8 


14.0 


2.25 


1.63- 


10^ 


40^ ^40 ^ g + 


1.321 


1.28 


1.31 


2.69- 


10^ 


^^K + e^^^ Ar + u 


1.513 




1.51 


3.33- 


103 



TABLE II. Masses, heat and neutrino production rates. M, H and L are in units of 
10^^ kg, 10^^ W and lO^^s"^ respectively. Hnr is defined as the difference between total heat flow 
and the radiogenic production. 



Model 


Chondritic 


BSE 


Fully Radiogenic 


M{U) 


0.65 


0.84 


1.7 


MiTh) 


2.5 


3.4 


6.5 


M{K) 


4.6 • 10"^ 


0.84 • 10"^ 


1.7- 10-4 


H{U) 


6.2 


7.9 


16.3 


H{Th) 


6.7 


8.6 


17.6 


H{K) 


16.4 


3.0 


6.1 


Hnr 


10.7 


20.5 





Lp{U) 


4.8 


6.2 


12.7 


UiTh) 


4.0 


5.2 


10.6 


MK) 


123 


22.5 


46.0 




15.2 


2.8 


5.7 



11 



TABLE III. Mass 


distribution and (ant i) neutrino fluxes. 


M and $ are in in units of 


10-*^^ kg, and 10® cm~^ s 


^ respectively. Uranium mass in 


the crust is 


fixed at the value estimated 


from rcf. [13]. (I>,y(rcac1 


tor) corresponds to the flux from 


a nuclear reactor with Pfj, = 2.8 GW at 


100 km. 








Uranium 


Model 


Chondritic 


BSE 


Fully Radiogenic 


M(crust) 


0.42 


0.42 


0.42 


M (mantle) 


0.23 


0.42 


1.29 


$p (crust) 


2.1 


2.1 


2.1 


$p (mantle) 


0.5 


1.0 


3.0 


*i7 (tot) 


2.6 


3.1 


5.1 



$,,(r('ac(<)r) {E < 3.2G McV) 0.1 



Thorium 



Model 


Chondritic 


BSE 


Fully Radiogenic 


M(crust) 


1.6 


1.6 


1.6 


M (mantle) 


0.89 


1.6 


4.9 


$p (crust) 


1.8 


1.8 


1.8 


$i7 (mantle) 


0.4 


0.8 


2.5 


(tot) 


2.2 


2.6 


4.3 



(reactor) {E < 2.25 MeV) 0.3 



Potassium 



Model 


Chondritic 


BSE 


Fully Radiogenic 


M(crust)/10^ 


0.42 


0.42 


0.42 


M(mantle)/10^ 


4.15 


0.42 


1.29 


(crust) 


7.7 


7.7 


7.7 


$p (mantle) 


34.9 


3.5 


10.9 


*p (tot) 


42.6 


11.2 


18.6 


$1, (crust) 


0.95 


0.95 


0.95 


$1, (mantle) 


4.33 


0.44 


1.04 


*v (tot) 


5.28 


1.39 


1.99 


(reactor) {E < 1.31 MeV) 




0.2 
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TABLE IV. Flux dependence on uranium abundance in the crust and distribution 
in the mantle. Fluxes arc calculated for Mi^^y^gfj^jyinjin^ (U) = 0.84 • IQi'^ kg. The total flux $(tot) 
is the sum of the contribution from the crust, upper mantle and transition zone+lower mantle. 
Same units as in Table III. 



Uranium 




high Ucc 


low ttcc 


depleted upper mantle 


M (crust) 


0.42 


0.21 


0.42 


M (upper mantle) 


0.06 


0.09 





M(transition+lower mantle) 


0.36 


0.54 


0.42 


$ (crust) 


2.1 


1.0 


2.1 


$ (upper mantle) 


0.21 


0.32 





$ (transition + lower mantle) 


0.79 


1.2 


0.92 


$(tot) 


3.1 


2.5 


3.0 



TABLE V. Flux dependence on the detector location. Fluxes in unit of 10^ cm ^ s ^ 
are calculated for Mcrust{U) = M^antle{U) = 0.42 • 10^^ kg, Th/U = 3.8 and K/U = W^. The last 
column is the average flux for the same values of Mcrust{U) and Mj^antleiU) (BSE model). 





Kamioka 


Gran Sasso 


average 


Uranium 


4.0 


4.6 


3.1 


Thorium 


3.4 


3.9 


2.6 


Potassium (u) 


14 


17 


11 


Potassium {u) 


1.8 


2.1 


1.4 
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